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Abstract 

It is known that some high-energy processes have a t-channel singularity in the 
physical region. In this paper we show that this singularity is regularized if one 
takes into account the finite sizes of the colliding beams, i.e. the realistic situation 
which takes place at high-energy colliders. On this way we obtain the finite cross 
section which is linear proportional to the transverse sizes of the colliding beams. 
As an application of the above result we study the production of the single W boson 
at ^'^ fi^ colliders in the reaction /i"//"*" — > eUgW^ . For the total energy ^/s ^ 95 
GeV and reasonable parameters of the muon beams the "beam-size " dependent 
part of the cross section for this reaction is of the order of the standard cross section 
of the fi~n~^ fi^i'^W^ process. 
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1 Introduction 



It is known since early 60 's that some high-energy processes can have a t-channel sin- 
gularity in the physical region This situation can occur if initial particles in a given 
reaction are unstable and the masses of the final particles are such that the real decay of 
the initial particles can take place. 

Recently in the paper 0] it was stressed that this problem turned out to be a practical 
issue for the reaction e~h'eW~^ . It is noted in that the standard calculation 

of this cross section leads to the infinite result. Indeed, if the invariant mass of the final 
eUe system is smaller than the muon mass m, the square of the momentum transfer in 
the t-channel can be both positive and negative depending on the scattering angles. 
In the region of small defined by the inequalities —A < q^ < A and A <^ the main 
contribution to the discussed cross section is given by the diagram with the exchange of 
the muonic neutrino in the t-channel (see fig.l). The t-channel propagator brings the 
factor to the matrix element M of this process 

Moc^^ (1) 
which results in the power-like singularity in the standard cross section 

daoc f \M\^dq'^(xB= f = oo . (2) 

-A -A ' 

It is, therefore, necessary to find physically reliable method to produce definite prediction 
for the measurable number of events. 

In the paper (see also it was suggested to regularize this divergence by taking 
into account the instability of muons in the initial state. If this is the case, then the muon 
mass is the complex quantity with the imaginary part proportional to the muon width 
r. It is possible then to solve the energy-momentum conservation constraints with the 
conclusion that the square of the momentum transfer in the t-channel acquires imaginary 
part proportional to the muon width. This observation leads to the replacement 

q^ ^ q^ — i'J, 7 ~ mT. 

The divergent integral in (^ is therefore regularized 

/r^-i-z^. (3) 



-A ' -A 

Consequently direct calculation of the scattering cross section becomes possible. 

We note in this respect that if such a regularization is performed, it is possible to 
estimate the typical time necessary for this reaction to occur. It is not difficult to see 
that this time is of the order of the moving muon life time, i.e. much more larger than 
the time necessary for muon beams to cross each other. Therefore for the realistic muon 
colliders the results of ref.0] are not applicable and a new consideration is necessary. 

We have found that for the realistic muon colliders there exist much more important 
effect connected with the finite beam sizes. In this paper we are going to show that 
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accounting for the finite sizes of the colliding beams gives a finite cross sections for the 
processes like 

fi^fi^ e~VeX (4) 

with the t-channel singularity in the physical region. This is the main resuh of our paper. 
As an example, we reconsider reaction fi~fj^~^ —>■ eUgW^ and show that the actual cross 
section of this process for the total energy ^/s ~ 100 GeV is approximately 1 fb for the 
typical transverse beam sizes a ~ 10"'^ cm. 

Let us mention here that the effect of the finite beam size at the high-energy colliders 
is well studied both experimentally and theoretically (for the review see ref. 0). For the 
first time this beam-size effect (BSE) was observed at the VEPP-4 collider (Novosibirsk) 
in 1980-81 during the study of the single bremsstrahlung in the electron-positron collisions 
This year the BSE was observed at HERA in the reaction ep — > epy [^. In both 
cases the number of observed photons was smaller than it was expected according to 
the standard calculations. The decreased number of photons was explained by the fact 
that impact parameters, which gave essential contribution to the standard cross section 
of these reactions, were larger by 2-3 orders of magnitude as compared to the transverse 
beam sizes. From the theoretical point of view, the BSE represents a remarkable example 
of the situation where traditional notion of the cross section and the standard formulae 
for the number of events are not valid. 

Let us briefly outline the principal points of our consideration. The standard rules 
for calculations in high-energy physics correspond to the collision of plane waves which 
represent initial states of colliding particles. Instead of this we consider collision of the 
wave packets which correspond to the particle beams in the initial state. As a consequence 
the standard rules are modified (see detailed discussion in the section 3). For instance, in 
the eq.(^) the following replacement takes place: 

\M\' ^ Mj,M}, . (5) 

Here Mfi and M/j/ are the matrix elements with identical final but different initial states. 
The difference in the momenta of the initial states is of the order of 1/a for the beam of 
the transverse size a. Due to this fact the singular propagators 1/g^ (which are the reason 
for the pathological behaviour of the cross section) are now taken in the different points 
of the momentum space for the amplitudes Mfi and M/j/ : 

1 1 7Tl 

q'' — X + ze + X + %e a 

Consequently the infinite factor B is replaced by: 

f dq^ f dq^ 

^ " ' |g2 + ^e|2 ^ ~ y (g2- A + ^e)(g2 + A-ze) 



-A 



Below we show that 



^ 1 

E ~ - ~ — 7 
A m 

and that the main contribution to this integral is given by the pole of the neutrino prop- 
agator + A — ze = 0, i.e. by the real neutrino. This means that the effect under 
discussion is connected with the decay fi~ — > e~z/eZ^^ and subsequent scattering of the real 
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muonic neutrinos on the /i"*" beam. As discussed in detail below these neutrinos have a 
broad energy spectrum (see fig. 2) 

dN^ = N^f{x)dx (8) 

{x is the fraction of the fi~ energy carried by neutrino). In accordance with the eq.(|^), the 
total number of these neutrinos (per one muon) is proportional to the typical transverse 
beam size a: 

K — (9) 

cr 

where r = l/F is the muon life time, cr = 660 m. As a consequence the "non-standard" 
contribution to the cross section is also proportional to the transverse beam size 

'-''non-stand = A^jy y fix)a^^^wixs)dx oca. (10) 

Let us mention here that in all the previous cases when the effects of the finite beam 
sizes have been studied, the results depend logarithmically on the beam size: 

cr oc In a . (11) 

In contrast to this, in our case ( t-channel singularity in the physical region ) the cross 
section (|I^) is linear proportional to the transverse sizes of the colliding beams. 

The subsequent part of the paper is organized as follows: in the next section, as a 
realistic example of the processes with the t-channel singularity in the physical region, 
we consider reaction (^ and present its cross section; in the section 3 we prove the basic 
formula used in the section 2; section 4 is devoted to the particular case — the reaction 
fj,~fj.~^ — > ePgW^] in the section 5 we present our conclusions; Appendix is devoted to 
the detailed consideration of a model process with a t-channel singularity in the physical 
region. 



2 General case 

In this section we show how the cross section of the reaction /x"//"'" — > e~VeX can be 
calculated. We introduce the following notations: s = {j>i + P2Y = ^E"^ is the square of 
the total energy in the center of mass frame, m and F = 1/r are muon mass and width 
respectively, p\ = p\ = m?, Ps is the 4-momentum of the final e~Ue system, y = p'^/m'^, 
q = Pi — Ps = {uj, q) is the momentum transfer in the t-channel and x = qp2/piP2 ^ uj/E. 
From simple kinematics it follows that 

9 x(l — X — y) n ^^^^ 

1 — X 1 — X 

where is the component of the momentum q which is transverse to the momenta of 
the initial muons. Note, that 

to > for y <l-x (13) 

and that for = 

|q±| = q± = m^Jx{l -x-y). (14) 
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For ?/ < 1 — X (or to > 0) we write the cross section in the form 

d(T = (icTstand + C^CTnon-stand (15) 

where by definition dagtsmd corresponds to the region < — and (icinon-stand corresponds 
to the region 

-m^<q^<tQ. (16) 

We will show that the main contribution to rfcTnon-stand is given by the region of a very 
small values of inside the region (|16D : 

- A < < A, m/a < A < (17) 

where a is the typical transverse beam size. In the region (|17D the main contribution comes 
from the diagram with the exchange of the muonic neutrino in the t-channel (fig.l). Since 
for such the exchanged neutrino is almost real, the corresponding matrix element can 
be considerably simplified. 

First, the numerator of the neutrino propagator q can be presented in the form 

Q ~ IIwA(q) MA(q)- 

A 

Here U\{q) is the Dirac spinor for the massless fermion with momentum q and helicity A. 
Second, in the region of interest the upper block of the diagram (see fig.l) corresponds to 
the decay /i~ ez/et'^. In this decay only left-handed z/^ is produced. Therefore in the 
above sum the only term with A = —1/2 is relevant. 

As a result, in the region (|l^ we present the matrix element M in the form 

M = -M^^eu^,^ M.^^x (18) 

q^ + ie 

where M^^ev^Vf, is the matrix element for the yU~ decay and M^^^x is the matrix element 
for the Pfiii'^ —>■ X process. In both of these subprocesses we take equal to zero. 



e 




Figure 1: The Feynman diagram for the reaction yU+ — eu^X which gives the leading 
contribution in the region of small 
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Using this equation we express the non-standard contribution through the muon decay 
width r and the cross section a^^^x of the z/^/x"*" — > X process as: 

1 dq^ 

C^O-non-stand = - X m dV — — ■ da^^^-^x- (19) 
vr \q^\^ 

For unpolarized muon beams we have 

r 

dT = -{l-y){l + 2y)dxdyd(p (20) 
vr 

where (p is the azimuthal angle of the vector q. Let us call the coefficient in front of 
da^n^x in eq. (|l9|) as the number of neutrinos: 

dNy = - X m dT = ^ x(l - y){l + 2y)dxdydip (21) 



As it is clear from the eq.(|T9D, the standard calculation of the cross section turns out 
to be impossible due to the power-like singularity, since the point = is within the 
physical region for y < 1 — x. 

The main result of our investigation of the BSE in the above process can be formulated 
as follows: 

Accounting for the BSE results in the following treatment of the divergent integral in 
the eq.l^: 

dq^ a 



B=It^--^ = - I (22) 
iy I y± mJx\l — x — y) 



The exact expression for the quantity a will be given below (see eqs.(|4^)-(^5D). We just 
mention here that it is proportional to the transverse sizes of the colliding beams. For 
the identical round Gaussian beams with the root-mean-square radii 

(^ix = (^iy = i = 1, 2 

this quantity is equal to 

a = \^a±. (23) 

Hence, the quantity i? ~ a/m as it is noted in (0). The contribution ( ]2^ ) comes from 
the region (|1^. It is not difficult to estimate that the contribution from the rest part of 
the region ( |16|) is smaller. Indeed, its relative value is of the order of 

i / ^ + (24) 
B J |g2|2 ^ BJ |g2|2 A ^ ' 

—m? A 

Using eqs. (P2|)-(P^ and integrating the number of neutrinos ( PT] ) over y (in the region 
< y < 1 — x) and over ip, we arrive to the following spectrum of neutrinos: 

dN^{x) _ Txa ^^^^ ^^^^ _ 24 —f^ ,22^ 16 2 
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where r is the hfe time of the muon at rest, cr = 660 m. The function f{x) is normahzed 
to 1: 

1 

f{x)dx = 1. 

The plot of this function is presented in the fig.2. The total number of neutrinos is equal 

After the spectrum of neutrinos is obtained, the non-standard cross section for the 
reaction /x"//"*" — >■ eue^ is given by the equation 

'^'^non-stand 

dNi,{x) da^f,^x{xs) = f{x) dx da^^^x{xs) . (27) 

Subsequent integration over x can be performed without further difficulties. For the 
particular case X = such calculation is performed in the section 4. 



3 Derivation of the basic formula 



Here we are going to prove our result presented in the previous section (see eq. (|2^) ). 

To begin with, let us note that the standard notion of the cross section is an approx- 
imation itself. As is well known, it corresponds to the plane waves approximation for 
initial and final particles. In the real experiments the particles are confined to the beams 
of a relatively small size and it is the collision of such beams that leads to the measurable 
number of events. 

In order to arrive to a more general formulae for the scattering processes, it is necessary 
to describe the collisions of the wave packets instead of the plane waves. In view of the 
fact that the movement of the particles inside a beam is quasiclassical, simple and efficient 
technique for taking into account the beam-size effects in the actual calculations has been 
developed (for the review see 0). 

Below we present some results from the refs. 0, which are essential for our dis- 
cussion. For simplicity, we neglect the energy and angular spread of the particles in the 
colliding beams. 

Let us remind that in the standard approach the number of events is the product 
of the cross section cr and the luminosity L: 

dN = da L, da oc |M|^ L = v J ni{r, t) n2(r, t)d^rdt (28) 

where v = |vi — V2I = 2 for the head-on collision of the ultra-relativistic beams. The 
quantities ni{r,t) are the particle densities of the beams. 

The transformation from the plane waves to the wave packets results in the following 
changes. The squared matrix element |Mp with the initial state in the form of the plane 
waves with the momenta pi and p2 transforms to the product of the matrix elements Mfi 
and Mjj, with the different initial states: 

da oc |Mp ^ da{K) oc MfiM}^, . (29) 

Here the initial state \i) is the direct product of the plane waves with the momenta 
ki = pi + |k and k2 = p2 — , while the initial state is the direct product of 
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the plane waves with the momenta k'^ = pi — |k; and ^ = p2 + . Instead of the 
luminosity L the number of events starts to depend on the quantity 

L{q) = V j Hi {r,t) n2{r + Q,t)d^rdt (30) 

through the following formula 

= /^f^"^'''^^^'")^^^)- ^^^^ 
The characteristic values of k are of the order of the inverse beam sizes, i.e. 

1 

K ~ — . 

a 

Usually this quantity is much smaller than the typical scale for the variation of the matrix 
element with respect to the initial momenta. In this case we can put k = in da{K, ) 
which results in the standard expression for the number of events (|28|) . Otherwise, one 
should analyze the complete formulae which take into account the effect of the finite beam 
sizes. 

In view of the discussion given in the previous section, this is indeed the situation 
which occurs in our case. Now we want to show how the finite result for the number of 
events can be obtained starting from the complete formula (l3l|) . 

Let us first define the "observable cross section" by the relation (compare with the 
eq-dH) 

, _dN _ 

d<J — -J— — ttdstand + "O'non-stand (32) 

where L is the standard luminosity. By writing the number of events in such a way, we 
push the BSE to the quantity cicxnon-stand- 



The study of the matrix element of the discussed process ( |18D suggests that the only 
quantity sensitive to the small variation of the initial momenta is the denominator of the 
neutrino propagator for small values of g^. Consequently, the transformations (^)-(|3T|) 
reduces to the following modification: 

1 11 

(33) 



\q^\^ t + ie t' — ie 
where 

q' = iPi-P3y, t={h-p,)\ t' = {k[-p,f. (34) 
Let us expand t and t' up to the terms linear in n. This gives 

t = q^ -\ t' = q^ + X (35) 

where 

X = kQ, Q = -p3 + — pi. (36) 
In the center of mass frame of the colliding muons 

Q ~ -P3± = q± = |q±| n, n = (cosv9,sinv9,0). (37) 
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As a result, the quantity B in the eq. (22) transforms to (compare with the eq. (|^)) 



r d^K(pQ L{q)_ r df^ 

J (27r)3 L _y (g2 - A + le) {q^ + X - le) ' ^ ' 

It has aheady been mentioned that the above equation must be understood in the sense 
that we isolate the part, which does not possess the limit k ^ 0, to the quantity B. All 
other quantities are smooth in the above limit and hence do not present any problem. 
To perform the g^-integration in (|38D we note [| that in the region (|lj 



X = KU^J {qiy - (1 - ^ Ao = Kn g|[. (39) 
This approximation is justified because 

k'l < A 



m? m? 



— < 1. 

Then we extend the region of integration over up to ±oo. The error introduced by this 



procedure is of the order of {■m/a)/A <^ 1 (cf. eg .(^41)). Now the integral over along 
the real axis can be replaced by the integral over the contour C which goes around the 
upper half plane of the complex variable q^: 

d(f f d(f TTi ^^^^ 



c 

Here 



D J D -Ao + ie 



D = {q^-Xo + te) (q^ + Xo-ie). 

The result is provided by the pole q^ + Xq — ie = in the upper half plane. 

Further integrations are simply performed with the help of the exponential represen- 
tation 



oo 



ia{-Xo+ie) ^4^^ 







-Ao + ie 

Subsequent integrations over k, and g become trivial and we obtain 



oo 





,0 



After simple redefinition aq ^^ = g we get for the factor B: 



This completes the proof of the eq. (|2^). 

^The more detailed discussion of this integration can be found in the Appendix. 
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At high energy colhders the distribution of particles inside the beams can be often 
considered as Gaussian. In this case L{gn) equals: 

r/ X T f 2/COS^(j9 Sin^09\1 , . r,\ /.ON 

L(gn) = Lexpj - g n = (cosy?, sin(/?, 0) (43) 

X y 

where 



2 2,2 2 2,2 



This results in the following expression for a: 

TT ay ^^^^ 

2 cos^ + a^. sin^ y? 

For the round (but not identical) beams with the root-mean-square radii 
we have 

TT 



2 V^?^ + ^2\- (45) 

It is interesting to note that the quantity a and the non-standard cross section are 
determined by the size of the largest beam. For the round and identical beams with 
ai±_ = CT2_L = a± the result of ( p3D can be obtained. 

4 The cross section of the eUeW^ reaction 

Let us analyze a special example taking X = . The cross section for the reaction 
z/^/i^ — > can be written in the form 



a^,^^ = Utt' ^^^ ^ ^""h ixs - M') (46) 
where r(iy /xz/^) = 0.22 GeV is the partial W decay width and M = 80.2 GeV is the 



W boson mass. Integration over the fraction of the neutrino energy x in the eq. (|27|) 
becomes trivial and we obtain: 

0^non-stand(/U"/U+ ^ CUeW^) = ^ aQ Xq /(Xq), Xq = — , (47) 

2CT s 

an = — = 20 nb. 

° M2 M 

For numerical estimates we take a = ^piio ^ (which corresponds to the case of the 
round identical Gaussian beams) with (see ref.0) 

= 10"^ cm. 

This non-standard cross section reaches the maximum of 0.76 fb for y/s = 93 GeV. For 
larger energies this cross section decreases as It is interesting to note that the 
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modest value of the non-standard cross section 0.76 fb is the result of the product of the 
very small number of neutrinos 

= = 4.2 ■ 10"^ 

2 cr 

and the huge value of the cross section for the z/^/U"*" transition averaged over 

neutrino spectrum 

{au^^w) = o-Q xq f{xo) = 1.8 ■ lO'^ fb, xq = 0.74. 

First, let us compare this non-standard piece with the standard contribution to the 
same cross section. We remind that by "standard" contribution we mean the cross section 
of the same reaction calculated by the standard rules excluding the region of the final 
phase space where > — m^. This contribution was calculated ^ with the help of the 
CompHEP package [|l^. The comparison of both contributions is shown in the fig.3. It 
is seen that the non-standard contribution dominates up to the energies -v/i ~ 105 GeV. 

Second, we compare our non-standard cross section with the cross section for the 
single W boson production in the reaction —>■ ^"v^W^ . The latter is completely 

standard process since it has no t-channel singularity. The reasonable estimate of its cross 
section can be quickly obtained with the help of the equivalent photon approximation. 
It gives the cross section 1 fb at ^/s ~ 95 GeV (where it almost coincides with our 
non-standard cross section). At higher energies the process fi^u^W~^ dominates 

as compared with the discussed process — > eVeW^ . 

5 Conclusions 

We have considered the processes with the t-channel singularity in the physical region. 
As the result of our study the following statement is proved: 

Accounting for the finite beam sizes eliminates the singularity and makes it possible to 
obtain the finite result for the observable cross section. 

In contrast to all previously known cases, the beam-size effect for the processes with 
the t-channel singularity in the physical region turns out to be linear proportional to the 
transverse size of the colliding beams. 

The effect appears to be connected with the real decay of the unstable particle in 
the initial state followed by the scattering of its decay products on the opposite beam 
1^. Evidently, the "interaction region " for such a scenario is only bounded if the finite 
sizes of the colliding beams are taken into account. In the traditional approach to the 
scattering process the initial states are the plane waves, which are not restricted to any 
finite volume. As a price, infinite cross section due to the t-channel singularity is obtained. 

The results presented in this paper correspond to the situation when the parameter 
a /{or) ~ Ni, <^ 1 or when the unstable particle life time is much larger than the interaction 
time ( for muon colliders this approximation is clearly perfect). If this parameter increases, 
stationary approximation becomes wrong at some point. In this case the effects of the 
initial particle instability and the beam size effects should be considered simultaneously. 

As a practical example, we consider the process of the single W boson production at 
fj^^fi^ colliders. We demonstrate that the standard cross section for a single W boson 
production and the non-standard one are of the same order of magnitude in the energy 
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region slightly above threshold of this reaction. For higher energies the relative size of the 
non-standard contribution rapidly decreases. 

We also note that precisely in the same way as above the processes like fj^^fi^ —>■ eu^W~^ 
and fi^^^ v^VeW^ can be considered. In this processes the t-channel singularity will 
be due to the electron neutrino and electron propagators respectively. In virtue of the 
discussion given above these processes will correspond to the collision of the real electron 
neutrino or electron produced in the /x~ decay with the opposite /x+ beam. 
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Appendix: The model example 

In this appendix a simple example of the reaction 2^2 with the t-channel singularity 
in the physical region is considered in detail. The same notations as in the main text of the 
paper are used. We work within the frame of the scalar field theory with the interaction 
Lagrangian 

Ci{x) = gA{x)B{x)C{x). 

The fields A, B and C represent the particles with the masses M, m and zero respectively. 
For definiteness we assume M > m. To avoid confusion, note also that the coupling g 
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has a dimension of mass. The particle A is unstable due to the A ^ B + C process, its 
decay width is of the order 

3 (48) 



IGttM 

Let us consider the scattering of the particles A and B : 

A{p,) + B{p2) ^ B{p;) + A{p^). (49) 

The matrix element of this process corresponds to the Feynman diagrams with the C- 
particle exchange in the t- and s-channels. For high energies, s S> M^, m^, the dominant 
contribution comes from the diagram with the C-particle exchange in the t-channel: 

M-. 



+ ie 



Simple kinematical relations provide the following limits for the square of the momen- 
tum transfer in the t-channel: 

(50) 

s 

Hence, the value of goes through zero. As was already indicated in the paper, in this 
case the cross section can not be calculated by means of the standard methods. It is easy 
to see, that = for the value of the transverse momentum 

|q±|=g°=V^= (51) 

V ^ 

To model the real situation discussed above, the following conditions have to be ful- 
filled: 

(i) The unstable particle life time must be much larger than the interaction time, i.e. 
(compare eqs. (|^) and (P^) 

MT < (52) 
a 

(a) The transverse beam size a must be large in comparison with the typical inverse 
transverse momentum ([51|), i.e. 

a > (53) 
v^o 



Taking into account (^81) and introducing small parameter 5 we express both of these 
constrains in the form 

T^T « 5 = ^ « 1- (54) 

We now apply the same technique, as in the main text of the paper. We start from the 
complete description of the scattering process in terms of the colliding beams (see eqs. (|29|)- 
(^)) and calculate the quantity (t(k). For this aim we expand t-channel propagators up 
to the terms of the order k, exactly as it has been done in the main text. Then the 
following formulae are valid: 

2 2 

Mfi = Mfi, = — (55) 

^ t + ie' ^ t' + ie ^ ' 
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where the quantities t and t' are given in (|3^)-(p7D. To calculate 



d(p J 327r2s2 ^ 

—s 

we have to calculate the integral 

^ = / D{q^) = it + ie){t' - ze) = (g^ _ ^q^ + te){q^ + Kq^ - te). 

Here the quantity 

nq± = fcn^to — (56) 

depends both on and y?. To evaluate J in the limit of the small momenta k, we rewrite 
it in the following way: 

f dq"^ ~f dq^ J dq^ 



Z}(g2) y D(g2) y D(g^ 

C — oo to 

where the contour C goes around the upper half planeQ. 

The integral over the contour C can be evaluated by using Cauchy's theorem. The 
result is provided by the pole in the upper half plane 

dq^ 2111 



c 



D{q^) R 



(58) 



where 

dD{q^ 



R = ^ at + Kq_L - ie = 0. (59) 

dq^ 

Taking into account the terms of the order 6 and neglecting the terms of the order 6^ , we 
obtain 

R = 2{-Kn^/t^ + ie). (60) 

As for the integrals both from to to the oo and from the — oo to — s, we note that the 
dependence of the integrand on k can be neglected and consequently 

' "1' (61) 



— ^ to 

As a result, we obtain 

dij{n) 9^ ( Tri 1 \ 

dip 327r2s2 /tny^ + ie to/ 
Using representation (^Tj) we get after trivial integration over k and Q 

do{AB BA) _ f na 1 \ 

d^ ~ \~^~To) 

where the quantity a is defined in (^). Generally speaking, a can depend on the angle 
(p (see (0)). Note that the "beam-size" dependent part dominates in the cross section 
(|63|) due to the inequality (|53|) . 



(62) 



(63) 



^Notc, that due to the eq.(p6|) the last integral m the eq.(p7|) also reqmres some definition. However 
for our purposes the dependence of this integral on k is not important. Hence the last integral in the 
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eq.(B^) should be understood as the value of this integral for 
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Figure 2: The normalized spectrum of neutrinos f{x) = {l/Ny) (dN^/dx) vs. x — the 
fraction of the muon energy carried by neutrino (see eq.(25)). 
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Figure 3: Standard (solid line) and non-standard (dashed line) contributions to the cross 
sections (fb) of the reaction fi~fi~^ ev^W^ in dependence on the total cms energy. The 
standard contribution is evaluated with the cut —(f> m^. 
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